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METHOD AND APPARATUS FOR AUTOMATED SIMULATION AND DESIGN OF 
CORNEAL REFRACTIVE PROCEDURES 

BACKGROUND OF THE INVENTION 

The invention relates to systems and techniques for 
5 mathematically modeling a human eye using topographical data 
measured from a human eye, using the mathematical model to 
simulate deformation of the eye by hypothetical incisions or 
excisions or ablations to arrive at an optimum surgical design 
by identifying the number, shape, location, length, and depth 

10 of the incisions or excisions or the shape, location, and 
quantity of corneal ablation ( s) . (It should be understood 
that hereinafter, including in the claims, the term 
"incision," which usually refers to a cut made by a scalpel, 
and the term "excision," which usually refers to a cut made by 

15 a laser beam, are considered to be interchangeable and to have 
the same meaning.) 

Modern corneal refractive surgery originated with the 
work of Dr. Svyatoslav Fyodorov of Moscow and Dr. Jose 
Barraquer of Bogota, Columbia. Subsequently, various surgical 

20 techniques have been developed to alter the curvature of the 
cornea to correct refractive errors. The various techniques 
include incisional keratotomy using diamond blades, excisional 
keratotomy using laser beams to photo-disrupt molecules and 
ablate tissue in a linear pattern, ablative keratectomy or 

25 photo-refractive keratectomy using laser beams to remove 
larger areas of corneal tissue, mechanical removal and 
reshaping of corneal tissue (keratomileusis) , and implantation 
of human or synthetic materials into the corneal stroma. All 
of the known procedures alter central corneal curvature by 

30 changing the structure of the cornea. All such refractive 
procedures are characterized by difficulty in predicting both 
the immediate arid long term results, because of errors in 
calculations of pre-surgical measurements, failure to 
precisely implement the planned surgical techniques, and 

35 biological variances which affect immediate and long term 
results . 
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The cornea traditionally has been treated as a 
spherocylindrical lens, assuming that the radius of each 
individual meridian from the corneal apex to the corneal 
periphery is uniform. Prior methodologies tend to use an 
5 approximation to the topographic information of the cornea to 
determine the refractive power of the eye. In one known 
procedure, circular mires (reflected light images from the 
cornea conventionally used to mathematically calculate corneal 
curvature) are reflected from the corneal surface, and the 

10 difference between a given point on the mire and an adjacent 
mire is measured. A semi-quantitative estimate of the surface 
curvature is obtained by comparing this measurement with the 
values obtained using spheres of various radii. Prior 
mathematical models use a variety of approximations such as a 

15 simplified form of the corneal surface (e.g., spherical) or 
assume a symmetrical cornea (leading to a quarter model or an 
axisymmetric model) or use simplified material properties 
(e.g., isotropic), or assume small deformations or 
displacements, or do not consider clinically obtained data in 

20 the construction of the mathematical model. 

Perhaps the closest prior art is indicated in the article 
"On the Computer-Aided and Optimal Design of Keratoref ractive 
Surgery," by Steven A. Velinsky and Michael R. Bryant, 
published in Volume 8, page 173 of "Refractive and Corneal 

25 Surgery," March/April 1992. This article describes a 
computer-aided surgical design methodology, proposing that it 
could be an effective surgical design aid for the refractive 
surgeon, wherein the surgeon could choose constraints on 
surgical parameters such as minimum optical zone size, maximum 

30 depth of cut, etc., measure the patient's corneal topography, 
refractive error and possible other ocular parameters, and 
then review the computed results. The article refers to 
several mathematical models described in the literature, and 
how such mathematical models might be helpful. However, the 

3 5 article fails to disclose any particular adequate mathematical 
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model of the cornea or any specific recommendation of surgical 
design that has been validated with clinical data. 

The prior radial keratotomy procedures frequently result 
in large amounts of undercorrection or overcorrection. Prior 
5 keratotomy procedures often are based on experiential use of 
nomograms indicating appropriate surgical designs for a 
particular patient based on age, sex, refractive error, and 
intraocular pressure. 

Finite element analysis is a known mathematically-based 

10 numerical tool that has been used to solve a variety of 
problems that are described by partial differential or 
integral equations. This technique has been used primarily in 
the area of solid mechanics, fluid mechanics, heat transfer, 
electromagnetics, acoustics, and biomechanics, including 

15 designing remedial techniques being developed for the human 
eye, to model internal structure and stresses in relation to 
various configurations of intraocular devices and corneal 
implants, as described in "Intraocular Lens Design With 
MSC/pal," by A.D. Franzone and V. M. Ghazarian in 1985 at the 

20 MSC/NASTRAN User's Conference in Pasadena, California, and in 
"Corneal Curvature Change Due to Structural Alternation by 
Radial Keratotomy," by Huang Bisarnsin, Schachar, and Black in 
Volume 110, pages 249-253, 1988 in the ASME Journal of 
Biomedical Engineering. Also see "Reduction of Corneal 

25 Astigmatism at Cataract Surgery," by Hall, Campion, Sorenson, 
and Monthofer, Volume 17, pages 407-414, July 1.991 in the 
Journal of Cataract Refractive Surgery. 

There clearly is an unmet need for an improved system for 
accurately predicting outcomes of hypothetical surgical 
30 procedures on the cornea to aid in the design of minimally 
invasive corneal surgery. There is an unmet need for a 
totally automated way of determining an optimal design of a 
surgical plan for incisional or excisional keratotomy or 
ablative keratectomy surgery to meet predetermined visual 
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objectives with minimum invasiveness and minimum optical 
distortion. Further, it would be desirable to provide a 
technique for designing a multi-focal cornea that is similar 
to a gradient bifocal for patients that have presbyopia. It 
5 would be desirable to have an accurate mathematical model of 
the cornea for use in developing new surgical procedures 
without experimenting on live corneas. 

SUMMARY OF THE INVENTION 

Accordingly, it is an object of the invention to provide 
10 a minimally invasive surgical procedure for corneal surgery 
for a human eye to achieve predetermined modified 
characteristics of that eye. 

It is another object of the invention to provide a system 
and method which result in improved predictability of outcomes 
15 of corneal surgery. 

It is another object of the invention to provide an 
improved method and apparatus for design of optical surgery 
that minimize invasiveness of the surgical procedure. 

It is another object of the invention to provide a method 

2 0 and apparatus for surgical design which result in reduction or 

elimination of postoperative irregular astigmatism. 

It is another object of the invention to provide an 
improved apparatus and method for surgical design which result 
in reduced multifocal imaging of the central cornea, thereby 
25 enhancing contrast sensitivity and improving vision under low 
light illumination conditions. 

It is another object of the invention to provide an 
improved finite element analysis model of the human eye, 
including back-calculation of values of material properties of 

3 0 the cornea and sclera, which incorporate the measured 
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topography of that eye and more accurately predict 
deformations of the cornea due to a hypothetical group of 
modeled incisions and/or excisions and/or ablation than has 
been achieved in the prior art. 

5 It is another object of the invention to provide a system 

and method for providing an optimal surgical design for a 
human eye to achieve desired optical characteristics thereof. 

It is another object of the invention to reduce the 

likelihood of postoperative complications in the eye 

10 including, but not limited to overcorrection or 
undercorrection of pre-existing refractive errors. 

It is another object of the invention to provide a 
"training tool" or "surgery simulator" for surgeons who need 
to gain experience with corneal refractive surgery. 

15 It is another object of the invention to provide a means 

for designing new surgical procedures without the need for 
experimentation on live human beings. 

Briefly described, and in accordance with one embodiment 
thereof, the invention provides a system for simulating 

2 0 deformation of a cornea as a result of corneal incisions, 
excisions and ablations in order to effectuate automated 
"surgical design" of a patient's eye in response to 
topographic, ultrasonic, and other measurements of the 
patient's eye. A finite element analysis model of the eye is 

25 constructed. Measured x,y,z coordinate data are interpolated 
and extrapolated to generate "nearest-fix" x,y,z, coordinates 
for the nodes of the finite element analysis mode. Measured 
thicknesses of the eye are assigned to each element of the 
finite element model. Pre-operative values of curvature of 

30 the cornea are computed. In one embodiment of the invention, 
Young's modulus is "back-computed" from measured values of 
corneal deformations at different pressure loads. An initial 
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estimated surgical plan, including a number of incisions, 
locations of incisions, incision orientations, incision depth 
and incision lengths is introduced into the shell finite 
element analysis model by introducing duplicate "nodes" and 
5 nonlinear springs along the initial hypothetical incisions. 
Or, ablations may be included in the estimated surgical plan 
introduced into the finite element analysis model by varying 
the thickness and/or material property constants of the 
elements in the ablated region. A geometrically and 

10 materially nonlinear finite element analysis then is performed 
by solving the equations representing the finite element 
analysis model in response to incremental increases in 
intraocular pressure until the final "equilibrium state" is 
reached. Postoperative curvatures of the cornea are computed 

15 and compared to pre-operative values and to vision objectives. 
If the vision objectives are not met, the surgical model is 
modified and the analysis is repeated. This procedure is 
continued until the vision objectives are met. In one 
embodiment, a boundary element analysis model is used instead 

20 of a finite element analysis model. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram illustrating the components 
used in the invention. 

Fig. 2 is a basic flow chart useful in describing the 
25 method of the invention. 

Fig. 3 is a block diagram of a subroutine executed in the 
course of executing block 35 of Fig. 2 to interpolate and 
extrapolate data in order to obtain the nodal coordinates of 
a finite element analysis model. 

30 Fig. 4 is a block diagram of another subroutine executed 

in the course of executing block 35 or Fig. 2 to "construct" 
the finite element analysis model. 
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Fig. 5 is a three-dimensional diagram of the finite 
element mesh used in accordance with the present invention. 

Fig. 6 is a partial side view illustrating both initial 
topography values of a portion of the cornea and final 
5 topography values resulting from simulated radial incisions 
and computed in accordance with the present invention. 

Fig. 7 is a diagram useful in explaining how incisions 
are included in the finite element analysis model of the 
present invention. 

10 Fig. 7A is a diagram useful in conjunction with Fig. 7 in 

explaining modeling of incisions. 

Fig. 8 is a diagram useful in explaining a technique for 
cubic spline interpolation and extrapolation to create 
"smoothed" three-dimensional data points from raw data 
15 provided by a keratoscope. 

Fig. 9 is a diagram useful in explaining automated back- 
calculation of the modulus of elasticity of the eye. 

Fig. 10 is a diagram useful in explaining optimization of 
the surgical plan according to block 41 of Fig. 2. 

20 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The present invention involves constructing a model of a 
human eye using a suitable three-dimensional finite element 
analysis (FEA) model which includes a mesh that generally 
corresponds to the shape of the human eye. The shape of the 
25 human eye is described by the finite element mesh obtained 
using data measured from a corneal mapping device and 
translated into the nodal points of the FEA model. The nodal 
points in a small region are connected to each other, to form 
a finite set of elements. The elements are connected to each 
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other by means of sharing common nodes. The thickness values 
at any particular region are obtained using ultrasonic 
measuring devices, and are appli.ed to the elements. The 
"loading" of the finite element mesh structure is represented 
5 by the intraocular pressure, and the resistance of the 
structure to such applied "loading" is measured by the 
stiffness of the structure, which is computed on the basis of 
its geometry, boundary conditions, and its material 
properties, namely Poisson's ratio v, and Young's modulus E. 

10 In the area of structural mechanics, finite element 

analysis formulations are usually based on the "principle of 
virtual work," which is equivalent to invoking the stationary 
conditions of the total potential energy, 11, given by 



2 v"" ^ v'" ^ s ^ '~ 



where 
(2) 

.t = B Z 

and 
(3) 



jl is the transpose of the strain vector, Z. is the transpose 
of nodal displacement vector, and Z^^ is the transpose of the 
nodal displacement vectors on the surface. and _Z^ are nodal 
displacement terms associated with nodal loads. In the above 
equations, the various symbols have the following meanings: 

represents the strain vector 
D represents the material matrix 
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Z represents the vector of nodal displacements 
represents the nodal body force vector 
represents the nodal surface traction vector 

dV represents differential volume 

dS represents differential surface area 

a represents the stress vector 

B represents a strain-displacement matrix 

V represents volume 

S represents surface area. 



10 The first term on the right hand side of the equation (1) 

is the strain energy of the structure, and the second and 
third terms represent the total work accomplished by the 
external forces and body forces. The strain energy is a 
function of the strains and stresses that are related to each 
- .15 other via the material matrix D. The material properties that 
contribute to the material matrix D include the modulus of 
elasticity (Young's modulus) and Poisson's ratio. In a 
uniaxial state of stress, Poisson's ratio is defined as: 



(4) 



20 where €iat is the normal strain in the lateral direction and 
eiong is the normal strain in the longitudinal direction. In 
a uniaxial state of stress, Young's modulus E is defined 
according to 

(5) 



where a„ is the normal stress and is the normal strain. 
The work accomplished is a function of the applied loads and 
surface tractions.. 
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Using an assumed displacement field, the minimization of 
the total potential energy n leads to the element equilibrium 
equations of the form 

(6) 



5 where the expression 

k = Im^ EEdV 

axn ^ 



is the element stiffness matrix, and z„«i is the vector of 
element nodal displacements is the vector of element nodal 
forces. Since the entire structure is assumed to be in 
equilibrium, the assembly of the element equations leads to 
10 the structural equilibrium equations of the form 

(7) 



Where K^xn is the structural stiffness matrix, Zn^i is the 
vector of nodal displacements and Rn^i is the vector of nodal 
forces. These algebraic equations are finally solved f or in 
15 a variety of ways depending on whether the structural behavior 
is linear or nonlinear. 

A commercially available finite element analysis program 
which effectively solves these equations after the appropriate 
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values and boundary conditions have been assigned to the 
various nodes and the appropriate material properties have 
been assigned to the various elements defined by the 
connectivity of the nodes is called ABAQUS, available from 
5 HKS, Inc. of Providence, Rhode Island. Creating the FEA model 
for purposes of the present invention simply involves 
inputting to the ABAQUS program the x,y,z coordinates for each 
node, inputting the forces that act on the nodes and/ or 
elements, assigning the boundary conditions to each node, 
10 defining the nodal connectivity that defines each element, and 
inputting the eye material properties and thickness or 
stiffness to each defined element along with other input data, 
such as whether the analysis is linear or non-linear, or the 
properties and definitions of the non-linear springs. 

It should be noted that there are two popular approaches 
to solving finite element analysis problems, one being the 
above-described approach of minimizing total potential energy 
(or, variational approach) , the other being a method of 
weighted residuals which operates on partial differential 
equations defining the problem. The first approach is 
generally recognized to be simpler, and is implemented by the 
above ABAQUS program, but the invention could be implemented 
by the second approach. 

Fig. 1 shows an apparatus used in conjunction with the 
25 present invention. An ultrasonic instrument 15, such as a DGH 
packymetor, model DGH-2 000 available from DGH Technology, Inc. 
of Frazier, Pennsylvania, is used to obtain the thickness and 
intraocular pressure of cornea llA. 

A corneal topographer 12, which can be a model TMS-1, 
30 manufactured by Computed Anatomy, 28 West 3 6th Street, New 
York, New York, is utilized to measure the surface topography 
of cornea llA. The resulting information is transferred by 
means of a floppy disk to a computer system 14. 
Alternatively, a digital data bus 13 could be provided to 
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transfer topography information from corneal topographer 12 to 
computer system 14. In a prototype of the invention, system 
14 includes an ordinary IBM PC type of computer 14A coupled by 
a RS232 cable 16 to a workstation computer 14 B, which can be 
5 a Sun SPARCStation 2, marketed by Sun Microsystems, of 
Mountain View, California. A printer 17 is connected by a 
cable to the printer port of computer 14 A. 

Thickness and intraocular pressure measurements made by 
ultrasonic instrument 15 are obtained. This data then is used 
10 in the generation of the finite element model. However, it is 
possible to have this data transferred digitally, either by 
means of a floppy disk or a communication link, to the 
personal computer 14 A. 

A conventional pressure loading device 19 is utilized to 
apply a precisely measurable force on a point of the sclera as 
far away as practical from the cornea, so that resulting 
changes on the elastic cornea as a result of the new loading 
can be measured. Then, in accordance with the present 
invention, the value of Young's modulus can be "back- 
calculated" in the manner subsequently described. 
Alternately, uniform pressure loading could.be achieved by 
applying a sealed pressure chamber to the eye and increasing 
the gas pressure therein. Such uniform loading may have the 
advantage of providing less "noise" error in the measurements. 
A suitable pressure loading device 19 could be an opthalmo 
dynamometer, commercially available from Bailliart, of 
Germany . 

To obtain an FEA model of the patient's eye, the measured 
topographical data is interpolated and extrapolated using the 
30 subsequently described cubic spline technique to provide a 
pre-established reduced number of nodal points of a finite 
element mesh, with nodal coordinates which are a "close fit" 
to the measured corneal surface. Values of the thickness of 
the cornea and sclera obtained from the data obtained from 



20 
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Ultrasonic instrument 5 are assigned to the various finite 
elements of the FEA model. The FEA mesh then defines a 
continuous surface that accurately represents the pre- 
operative surface of the cornea, including any astigmatism 
5 that may be present. 



The curvatures of the surface then are computed at each 
node of the finite element analysis model. Surfaces of 
revolution are generated by revolving a plane curve, called 
the meridian, about an axis not necessarily intersecting the 

10 meridian. The meridian (defined by a radial line such as 21 
in Fig. 5) is one of the principal sections and its curvature 
at any point is one of the principal curvatures k,. (The 
principal curvatures are defined as the maximum and the 
minimum curvatures at a point on a surface.) The other 

15 orthogonal principal section is obtained by the intersection 
of the surface with a plane that is at right angles to the 
plane of the meridian and that also contains the normal. The 
second principal section has curvature kz. If the equation of 
the meridian is written as r=f(z), then 



20 (8) 



[l+(r')^] ■ 



and 
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r' and r" being the first and second derivatives of r, 
respectively. 

Fig. 2 is a flowchart useful in explaining the basic 
steps involved in use of the system shown in Fig. 1 to produce 
5 an optimum design for guiding surgery of a patient's eye. In 
block 31 of Fig. 2, the physician determines the "vision 
objectives" for the eye. The vision objectives can be 
specified as desired curvatures at particular locations of the 
cornea when the eye is in an accommodatively relaxed condition 

10 after completion of the surgery. These desired curvatures of 
refractive change may be determined at nodal points on the 
cornea by use of a spatially resolved ref ractometer . The 
vision objectives are selected to maximize the number of light 
rays that the eye focuses on the fovea for a given functional 

15 distance. 

In block 32, the physician provides initial estimates of 
the number of incisions required, their locations, the 
orientations of various incisions, the incision lengths, and 

the incision depths needed in order to accomplish the vision 
20 objectives of block 31. As indicated in block 33, corneal 
topographer 12 is used to obtain a topographic map of a 
portion of the eye. The TMS-1 corneal topographer mentioned 
above is capable of providing an x,y,z coordinate "map" that 
covers most of the cornea, producing a data file from which 
25 the x,y,z coordinates of approximately 7000 points can be 
obtained . 

As indicated in block 34 of Fig. 2, the ultrasonic 
instrument 15 is used to provide measurements of the thickness 
of the cornea and the intraocular pressure. In the prototype 
3 0 system presently being implemented, typical values of 
Poisson's ratio and Young's modulus are used. Presently, 
Poisson's ratio values of 0.49 are used for both the cornea 
and sclera. Presently, values of Young's modulus equal to 2 
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X 10^ dynes per square millimeter are used for the cornea and 
5 X 10^ dynes per square millimeter for the sclera. 

Preferably, Young's modulus is "back-calculated" on the 
basis of corneal topographical changes measured by using the 
5 corneal topographer 12 after varying a known force applied by 
pressure loading device 19 (Fig. 1) to the eye. The main 
objective of the back-calculation procedure is to determine as 
accurately as possible the modulus of elasticity for the 
cornea and the sclera, because it also is recognized that 

10 these values vary from patient to patient, and because it also 
is recognized that the modulus of elasticity is one of the 
most crucial parameters that influences the finite element 
analysis predictions. To describe the basic technique of the 
back-calculation procedure, refer to Fig. 9, which shows three 

15 assumed states, namely State 0 in which the cornea is relaxed. 
State 1 in which pressure loading device 19 applies point load 
PI to the sclera of eye 11, and State 2 in which pressure 
loading device 19 applies point load P2 to the sclera. PO is 
the intraocular pressure which is uniformly applied to the 

20 inner surface of the cornea and sclera. The values of the 
moduli of elasticity for the cornea and the sclera, 
respectively, are adjusted such that the z coordinates at 
selected nodes are close to the values actually measured for 
State 1 and State 2 by TMS-1 corneal topographer 12 with the 

25 two values of point loads PI and P2 actually applied, 
respectively. 

Let Z'l be the observed z coordinate at a particular node 
i for a state j obtained using the above mentioned TMS-1 
system. Let Z~ be the computed z coordinate at a particular 
3 0 node i for the state j using the finite element analysis 
according to the present invention. The back-calculation 
problem then is to find the value of {E<,,EJ to minimize the 
value of the expression 
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with the conditions 



and 



and where {E<.,Es} is the vector of design variables, f(Ec,Es) 
5 is the objective function, E,. is the modulus of elasticity of 
the cornea, E^ is the modulus of elasticity of the sclera, n 
is the number of points at which z displacements are to be 
computed, and the two inequality constraints represent the 
lower (L) and upper (U) bounds on the two design parameters. 

10 It should be appreciated that such a problem formulation falls 
under the category of a non-linear programming (NLP) problem, 
and can be solved using various non-linear programming 
techniques such as the "method of feasible directions", or 
using a constrained least-squares technique. A commercially 

15 available program for solving such non-linear programming 
problems is the DOT (Design Optimization Tools) program, 
available from VMA, Inc. of Santa Barbara, California. 
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As indicated in block 35, an FEA model is "constructed" 
for the patient's eye by interpolating between the various 
7000 x,y,z coordinates of the corneal map produced by TMS-1 
corneal topographer 12 to provide a smaller number of 
5 representative "smoothed" x,y,z values to be assigned to the 
various nodes of the FEA mesh shown in Fig. 5. 

Fig. 5 shows one quadrant of the FEA mesh, the other 
three quadrants being substantially identical except for the 
nodal values assigned to the nodes thereof. The FEA mesh 
shown in Fig. 5 includes a plurality of equi-angularly spaced 
radial lines 21, each extending from a cornea center or apex 
27 of cornea section 24 to the bottom of sclera section 23. 
In the FEA mesh actually used in a prototype of the invention 
under development, there are 32 such radial lines 21 and also 
30 generally equally spaced circumferential lines 22. Each 
area such as 25 that is bounded by two adjacent radial lines 
and two adjacent circumferential lines 22 is an "element" of 
the FEA model. Each typical element 25 includes eight 
assigned "nodes", such as nodes 26-1... 26-8. The four corners 
of a typical element such as 25 share corner nodes 26-1,3,5,7 
with adjacent elements, and also share "midpoint" nodes 26- 
2,4,6, and 8 with corresponding midpoint nodes of adjacent 
elements. The nodes and the connectivity thereof which define 
the elements of the FEA mesh thus are illustrated in Fig. 5. 

25 The values assigned to each node include its 

interpolated/extrapolated x,y,z coordinates and dts boundary 
conditions, which are whether the node can or cannot undergo 
x,y,z displacements and rotations. The values assigned to 
each element in the FEA model include the thickness of the 

30 element. Young's modulus or modulus of elasticity, the shear 
modulus, and Poisson's ratio in the orthotropic directions, 
namely the xy, xz, and zx directions. Any external "loading" 
forces at each node also are assigned to that node. The 
orthotropic values of Poisson's ratio presently uses are v^= 



15 
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0.0025, Vx2= 0.0025, and v^^= 0.49. The value of shear modulus 
used is Gy^= 6.71 x 10^ dynes per square millimeter. 

The objective of the tasks in block 35 is three-fold. 
First, the total number of nodes, and thus the elements 
5 generated from them, should be a variable, so that the mesh 
sensitivity of the results can be studied while the operator 
is given the chance to use a coarse mesh for preliminary 
studies. Second, the nodal points generated should be 
compatible with the choice of element required. For example, 

10 eight-node shell elements are used in the present approach. 
However, the proposed system is able to generate any type of 
element required, such as a 27-node hexahedral three- 
dimensional element, a 6-node triangular shell element, or a 
9-node shell element. Third, the nodes generated must be able 

15 to provide sufficient mesh refinement or density to achieve 
the needed resolution. A refined mesh in the regions of 
primary interest such as the optical zone is important since 
it can capture the stresses and the variations of the 
displacements, and thus, the curvatures. The mesh refinement 

20 parameter is chosen by the operator as one of the variables to 
study for the regions of primary interest, such as the optical 
zone (i.e., the portion of the cornea central to the radial 
incisions) , while the other regions such as the sclera are 
still incorporated in the model. 

25 Some of the steps performed by computer 14 in accordance 

with block 35 of Fig. 2 are shown more specifically in Figs. 
3 and 4. As indicated in block 45 of Fig. 3, computer 14A 
reads the ASCII data files containing the above-mentioned 7000 
coordinates of the corneal map produced by the TMS-1 corneal 

30 topographer 12. Due to the nature of the data collection, it 
is possible that some "noise" exists in the original data. 
The origin of the noise might be attributed to the inability 
of corneal topographer 12 to provide an exact determination of 
the coordinates, or the lack of existence of the coordinate 

35 value at an expected site. As indicated in block 46, a simple 



wo 94/18636 



PCT/US94/00908 



19 

program scans the original data for elimination of such data 
points. 

As indicated in block 47 of Fig. 3, the polar coordinate 
data supplied by the TMS-1 corneal topographer is converted 
5 into the above-mentioned 7000 x,y,z coordinates. Points which 
lie along the radial lines 21 of the FEA mesh shown in Fig. 5 
are selected for use in the interpolation/ extrapolation 
process described below. 

As indicated in block 48, the cubic spline interpolation 
10 and extrapolation procedure (described later with reference to 
the diagram of Fig. 8) is utilized to compute the intermediate 
x,y,z coordinates for each node of the FEA mesh lying on the 
pre-defined radial line 21 (Fig. 5) . Then, as indicated in 
block 49 of Fig. 3, the program creates a final set of x,y,2 
15 coordinates for nodes that lie on the circumferential lines 22 
of the FEA mesh (Fig. 5) using the cubic spline 
interpolation/extrapolation method. This step is necessary 
since data points that have the same radial coordinate do not 
necessarily have the same height or z value. 

20 At this stage, as indicated in block 50 of Fig. 3, the 

initial or pre-surgery diopter values at the final setup 
points are computed. Once the radial lines 21 are 
"generated", a series of nodes are selected at a specific 
height and used to obtain the circumferential nodes, i.e., the 

25 nodes which are on circumferential lines 22 "between" the 
radial lines 21. The x,y, z' coordinates and diopter values of 
curvature at each node of the model then are output. 

As indicated in block 52 the coordinate data files 
produced by block 51 of Fig. 3 are read. In block 53 the 
30 finite element mesh options/data are read. This includes the 
orthotropic material properties of the cornea and sclera, the 
nonlinear load-elongation curve data used by the "spring" 
elements which represent incision depth as subsequently 



wo 94/18636 



PCT/US94/00908 



20 

described with reference to Figs. 7 and 7A, the loading 
information (i.e. , the intraocular pressure) , and the boundary 
conditions (i.e., the connections of the bottom nodes of the 
sclera to a stationary reference) . Then, the program reads 
5 the incision data, as indicated in block 54, and goes to block 
55 in which the FEA model is "created", i.e., duplicate nodes, 
element connectivity and load-elongation data for the spring 
elements are created. 

Finally, in block 56, the data files re'iquired for 
10 carrying out a geometric and materially nonlinear finite 
element analysis are created and output. In the present 
embodiment of the invention, the above-mentioned ABAQUS 
program is used as the finite element analysis program and is 
executed on the workstation computer 14 B. 

15 Returning to Fig. 2, in block 36, pre-operative 

curvatures are completed in diopters at each node of the FEA 
model . 

Then, in block 37, an initial (estimated) number of 
incisions are "constructed" in the FEA model using the 

20 information established in block 32. Figs. 7 and 7A 
illustrate how each such incision is modeled in accordance 
with the present invention. In Fig. 7, the FEA mesh llA 
includes an incision 61 modeled along a radial line 58 of the 
FEA mesh, in which numerals 61-1,2,3,4,5 represent all nodes 

25 of the FEA mesh from one end of the modeled incision 61 to the 
other. Additional "duplicate" nodes 62-2,3,4 are added to the 
FEA model, immediately adjacent to nodes 61-2,3,4, 
respectively. It should be noted that the "connectivity" of 
the nodes of the FEA mesh are omitted between corresponding 

30 duplicate nodes 62-1,2,3,4 and nodes 61-1,2,3,4, respectively. 

These duplicate nodes are connected with spring elements 
63-2,3,4 which have nonlinear load-deflection curves. The 
nature of the curves is a function of the depth of the 
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incision and the material properties of the tissue through 
which the incision is made. The depth of the modeled incision 
61 is represented by equations corresponding to the nonlinear 
elastic spring elements such as 53-1,2,3,4 in Fig. 7A, which 
are the only direct interconnections between the FEA mesh 
nodes and the immediately adjacent duplicate nodes. When an 
FEA program is executed, the effect of the intraocular 
pressure is to cause the incision 61 to "open" or spread an 
amount determined by the elastic spring constants assigned to 
the spring elements 63-1,2,3 and 4. 

The above-mentioned ABAQUS FEA program, when executed as 
indicated in block 37 of Fig. 2, computes the displacements at 
each node of the FEA model in response to the intraocular 
pressure . 

The computed nodal x,y,z displacements are added to the 
corresponding pre-operative x,y,z values for each node, and 
the results are stored in a data file. If desired, the 
results can be displayed in, for example, the form illustrated 
in Fig. 6. Post-operative curvatures (computed in diopters) 
then are computed for each node based , on the new nodal 
locations. 

In Fig. 6, which shows a computer printout produced by 
the system of Fig. 1, the measured pre-operative configuration 
of the eye surface is indicated by radial lines 82, and the 
25 computed post-operative configuration is indicated by radial 
lines 84. Numerals 21 generally indicate radial lines of the 
FEA model, as in Fig. 5. More specifically, numerals 21A-1 
and 21A-2 designate radial lines of the pre-operative surface 
represented in the FEA model, and numerals 2lB-l and 21B-2 
30 represent radial lines of the "computed" post-operative 
surface in the FEA model. Numerals 61A designate proposed 
radial incisions in the "measured" pre-operative surface, and 
numerals 61B designate the same incisions in the "computed" 
post-operative surface. (The individual circumferential lines 
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of the computer printout of Fig. 6 are difficult to identify, 
but this does not prevent accurate interpretation of the 
effect of the proposed incisions on the curvature of the 
cornea.) Numerals 86 indicate the limbus. 

5 As indicated in block 38, the post-operative curvatures 

then are compared with the pre-operative curvatures and the 
vision objectives established according to block 31 to 
determine whether the initial estimated surgical plan 
accomplished the vision objectives. 

10 Then, as indicated in block 39, computer 14 determines if 

the vision objectives are met. If the determination of 
decision block 39 is affirmative, the surgical design is 
complete, as indicated in label 40. Otherwise, however, the 
program executed by computer 14 goes to block 41, and an 

15 optimization technique, subsequently described with reference 
to Fig. 10, is utilized to modify the number of incisions and 
their locations, orientations, lengths, and depths. The 
process then returns to block 37 and repeats until an 
affirmative determination is reached in decision block 39. 

20 The technique for modifying and optimizing the surgical 

design according to block 41 can be understood with reference 
to Fig. 10. As indicated above, the vision objectives are to 
obtain prescribed curvature values at specific FEA nodal 
locations i on the cornea. As an example, assume the surgical 

25 plan includes the locations of the incisions and includes the 
location and depth of each incision. The surgical 
optimization problem of block 41 then can be defined to be the 
problem of determining a value of {aj,lj,dj) that minimizes the 
value of the expression 
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where {aj,lj,dj} is the vector of design variables, f{aj,lj,dj} 
is the objective function, aj is the starting radial distance 
from the center of the finite element model, as shown in Fig. 
10, Ij is the length of the incision, dj is the depth of the 
incision, j is the incision number, ri is the computed 
curvature value based on the results from the finite element 
analysis, r^ is the observed curvature, n is the number of 
points at which curvature computations are to be carried out, 
and the three above inequality constraints represent the lower 
(L) and upper (U) bounds on the three design parameters. 

With this information, it is possible to include any 
parameter that influences the finite element model as a 
potential design variable, and any response or parameters 
related to the response computed by the finite element 
analysis as appearing in the objective functions or 
constraints. For example, the shape of the incision, the 
number of incisions, or the thickness of ablative tissue can 
be design variables. 
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The foregoing problem formulation falls under the 
category of nonlinear programming problem. Those skilled in 
the art can readily solve such problems using nonlinear 
programming techniques utilizing commercially available 
5 nonlinear programming software, such as the previously 
mentioned DOT program. 

It should be understood that computer 14 includes not 
only IBM PC-type computer 14A (which can be a 86386 or 86486 
based computer) , but also a workstation computer 14B, merely 
10 because the ABAQUS program cannot be executed on a PC-type 
computer. If, in the future, a suitable FEA analysis program 
becomes available, all of the computation can be performed by 
an IBM PC-type computer 14A. 

Fig. 8 is useful in illustrating the application of cubic 

15 spline techniques referred to in block 48 of Fig. 2 to 
interpolate/extrapolate data from the nodal points of the FEA 
mesh from the data points obtained from TMS-l corneal 
topographer 12. In Fig. 8, numeral 71 designates the z axis 
or a center line of the cornea passing through its apex, 

20 numeral 72 designates a radial line along which data points 
obtained from corneal topographer 12 lie, and numeral 73 
designates various such measured data points. The extent of 
the cornea is indicated by arrow 78, and the extent of the 
sclera is indicated by arrow 79. The extent of the "optical 

25 zone" is indicated by arrow 80. As indicated above, the TMS-l 
corneal topographer provides 7000 such data points 73. The 
first step of the cubic spline process takes such data points, 
as indicated by arrow 74, and "fits" each segment of radial 
line 72 between adjacent corneal topographer data points 73 to 

30 the equation 



(12) 



z = ax^ + bx^ + cx + d. 
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where z is a distance along center line 71, and x is distance 
in the horizontal direction from line 71 toward the base of 
the sclera. Equation 12 then is used to compute values of z 
for each value of x corresponding to a node of the FEA mesh 
5 (shown in Fig. 5) to obtain values z for each of the nodes of 
the FEA mesh along each radial line 21, as indicated by arrow 
76 in Fig. 8. Values for the "midpoint" nodes such as 26-2 
and 26-6 of Fig. 5 are obtained by interpolating adjacent 
nodal values of z on the same circumferential line 22. Most 
10 texts on numerical analysis disclose details on how to use the 
cubic spline technique, and various commercially available 
programs, such as IMSL, available from IMSL, Inc. of Houston, 
Texas, can be used. 

A fixed boundary condition for the base of the sclera can 
15 be assigned. It has been found that the nature of the 
boundary conditions at the base of the sclera has only a small 
effect on the results of the finite element analysis of the 
cornea . 

The above-described model was utilized to compute the 
20 strains and nodal deflections in a particular patient's eye 
based on measured topographical data extending outward 
approximately 8 millimeters from the center of a patient's 
eye. The measured data was extrapolated outward another 8 
millimeters to approximate the topography of the remaining 
25 cornea. Computed results for six eyes of four different 
patients using the above FEA model agreed with clinically 
measured results to within one to ten percent. It is expected 
that when the back-calculation of Young's modulus is 
incorporated into the process, agreement will be improved, 
30 especially in the optical zone and close to incisions. 

The above-described FEA model cari- be used to pre- 
operatively design incisions, excisions or ablations of the 
cornea, resulting in great predictability of surgical outcome 
and thereby allowing minimum invasiveness to achieve the 
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desired result with the least amount of surgical trauma to the 
cornea. Fewer operative and post-operative visits by the 
patient to the surgery clinic are likely. Advantages of the 
improved surgical designs which result from the above- 
5 described invention include reduced multi-focal imaging of the 
central cornea, thereby enhancing contrast sensitivity and 
improving vision under low light illumination conditions. 
Reduction or elimination of post-operative irregular 
astigmatism is another benefit. Yet another benefit is 

10 minimization of side effects such as glare and fluctuation of 
vision associated with traditional incisional keratotomy. The 
described mathematical model will have other uses, such as 
allowing design of a bifocal corneal curvature to allow both 
near and distance vision for patients in the presbyopic stage 

15 of their lives. The model of the present invention also will 
allow development of new surgical techniques for correcting 
nearsightedness, farsightedness and astigmatism as a viable 
alternative to experimenting on live human corneas. 

While the invention has been described with reference to 

20 several particular embodiments thereof, those skilled in the 
art will be able to make the various modifications to the 
described embodiments of the invention without departing from 
the true spirit and scope of the invention. It is intended 
that all combinations of elements and steps which perform 

25 substantially the same function in substantially the same way 
to achieve the same result are within the scope of the 
invention . 

For example, keratoscopes or other cornea measurement 
devices than the TMS-l device can be used. Non-radial 

30 incisions, such as T-shaped incisions for correcting 
astigmatism, can be readily modeled. The stresses at various 
locations can be displayed, for example, by graded colors. 
Many variations of the finite element model are possible. In 
the two-dimensional shell finite element analysis model 

35 described above, the use of the nonlinear springs to model 
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depths of incisions could be avoided by modeling elements 
around the proposed incision to have reduced thickness and/or 
different material properties, so that the incision region has 
reduced stiffness, and the computed deformations are 
5 essentially the same as if the nonlinear springs were to be 
used. For example, it is possible to use three-dimensional 
finite elements in lieu of the two-dimensional shell finite 
elements with assigned thickness parameters, and model the 
incisions directly, without having to use the nonlinear spring 

10 elements. Mathematical models other than a finite element 
analysis model can be used. For example, a boundary element 
analysis model could be used. As those skilled in the art 
know, the basic steps in the boundary element methods are very 
similar to those in the finite element methods. However, 

15 there are some basic differences. First, only the boundary is 
discretized, that is, the elements are "created" only on the 
boundary of the model, whereas in finite element analysis 
models the elements are "created" throughout the domain of the 
model. Second, the fundamental solution is used which 

20 satisfies the governing differential equation exactly. A 
fundamental solution is a function that satisfies the 
differential equation with zero right hand side (i.e., with 
body force set to zero) at every point of an infinite domain 
except at one point known as the source or load point at which 

25 the right hand side of the equation is infinite. Third, the 
solution in the interior of the model can be obtained 
selectively once the approximate solution on the boundary is 
computed. Although constant intraocular pressure has been 
assumed, non-constant intraocular pressure could be 

3 0 incorporated into the described technique. Although post- 
operative swelling has been assumed to not effect the eventual 
curvatures of the cornea, healing of the incision does effect 
the eventual curvature. The finite element analysis model can 
be adapted to model such healing effects and predict the final 

35 curvatures, stresses, etc. 
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Additionally, p-finite elements, Raleigh-Ritz , mixed 
formulations, Reissner's Principal, all can be used to 
generate the finite element equations. These equations, then, 
can be used in the modeling method of the present invention. 
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WHAT IS CLAIMED IS : 



- 1 1. A method of simulating deformation of a cornea of an 

2 eye in response to a corneal incision, comprising the steps 

. 3 of: 

4 (a) measuring the topography of a portion of the 

5 eye using a topography measuring device to produce coordinates 

6 for a number of data points of a surface of the eye; 

7 (b) storing in a computer system a finite element 

8 analysis model of the eye, the finite element analysis model 

9 including a number of nodes, the connectivities of which 

10 define a plurality of elements; 

11 (c) operating the computer system to interpolate 

12 between and extrapolate beyond the data points to obtain a 

13 reduced number of x,y,z coordinates that correspond to nodes 

14 of the finite element analysis model, respectively, and 

15 assigning the x,y,z coordinates to the various nodes, 

16 respectively; 

17 (d) determining a value representing intraocular 

18 pressure in the eye and assigning that value to each element; 

19 (e) representing a first incision in the finite 

20 element analysis model by means of a plurality of duplicate 

21 nodes adjacent to a group of incision-bounding nodes of the 

22 finite element analysis model located along the first incision 

23 and representing the depth of the first incision by means of 

24 a plurality of nonlinear spring elements each connecting an 

25 incision-bounding node to an adjacent duplicate node, 

26 respectively; 

27 (f) using the finite element analysis model to 

28 compute, at each of the nodes, new values of the x,y,z 

29 coordinates resulting from the first incision; and 

'30 (g) displaying the nodes having the computed x,y,z 

31 coordinates to show the simulated deformation of the cornea. 



1 2. The method of Claim 1 including establishing a vision 

2 objective for the eye and storing vision objective information 
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3 in the computer system, wherein step (f) includes comparing 

4 the simulated deformation of the cornea with the vision 

5 objective to determine if the first incision results in the 

6 vision objective being met. 

1 3. The method of Claim 2 including, if the vision 

2 objective is not met, modifying the first incision and/or 

3 adding a second incision in the finite element analysis model 

4 similar to the first incision, and repeating step (f) to 

5 determine if the vision objective is met. 

1 4. The method of Claim 3 wherein step (c) includes 

2 executing a computer program containing the finite element 

3 analysis model so as to minimize the total potential energy of 

4 the surface of the eye represented in the finite element 

5 analysis model. 

1 5. The method of Claim 4 including measuring the 

2 thickness of various points of the cornea and/or sclera and 

3 assigning values of the measured thicknesses to each element 

4 of the finite element analysis model, respectively, before 

5 step (f ) . 

1 6. The method of Claim 4 including modeling an ablation 

2 of the cornea in the finite element analysis model by 

3 assigning a plurality of values of thickness corresponding to 

4 an ablated portion of the cornea to elements of the finite 

5 element analysis model in the ablated portion of the cornea, 

6 respectively. 
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1 7. The method of Claim 4 wherein step (e) includes 

2 representing a plurality of radial keratotomy incisions, 

3 including the first and second incisions, in the finite 

4 element analysis model. 

1 8. The method of Claim 4 wherein the first incision is 

2 a non-radial incision. 

1 9. The method of Claim 4 including assigning values of 

2 material constants of the eye, including Poisson's ratio, 

3 modulus of elasticity, and shear modulus, to each element of 

4 the finite element analysis model. 

1 10. The method of Claim 7 wherein the modifying includes 

2 executing a nonlinear programming computer program to 

3 determine how much to modify the number of incisions, the 

4 orientations of the incisions, and the lengths and depths of 

5 the various incisions. 

1 11. The method of Claim 2 wherein establishing the 

2 vision objective includes providing an initial set of surface 

3 curvatures for the cornea, the method including computing 

4 simulated post-operative curvatures from the new values of 

5 X/y>z coordinates computed in step (f) and comparing the 

6 simulated post-operative curvatures with the surface 

7 curvatures of the initial set to determine if the vision 

8 objective is met. 

^ 1 12. The method of Claim 2 wherein each element of the 

2 finite element analysis model is an eight-node element, and 

3 wherein a boundary condition of the finite element analysis 
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4 model is that a base portion of the finite element analysis 

5 model is stationary. 



1 13 . The method of Claim 3 including assigning 

2 substantially different measured values of modulus of 

3 elasticity to elements of cornea portions and sclera portions 

4 of the finite element analysis model. 



1 14. The method of Claim 1 wherein step (c) includes 

2 executing a cubic spline computer program to obtain the 

3 reduced number of x,y,z coordinates according to an equation 

4 z = ax^ + bx^ + cx + d which has been fit to the measured data 

5 points of step (a) , x being distance from an apex axis of the 

6 eye. 



1 15. The method of Claim 3 including selecting a vision 

2 objective which produces a simulated multi-focal configuration 

3 of the cornea. 



1 16. A method of simulating deformation of a cornea of an 

2 eye as a result of a corneal ablation, comprising the steps 

3 of: 

4 (a) measuring the topography of a portion of the 

5 eye using a topography measuring device to produce coordinates 

6 for a number of data points of a surface of the eye; 

7 (b) storing in a computer system a finite element 

8 analysis model of the eye, the finite element analysis model 

9 including a predetermined number of nodes, the connectivities 

10 of which define a plurality of elements; 

11 (c) operating the computer system to interpolate 

12 between and extrapolate beyond the data points to obtain a 

13 reduced number of x,y,z coordinates that correspond to nodes 



wo 94/18636 



PCT/US94/00908 



3 3 

14 of the finite element analysis model, respectively, and 

15 assigning the x, y, z coordinates to the various nodes, 
. 16 respectively; 

17 (d) determining a value representing intraocular 

18 pressure in the eye and assigning that value to each element; 

19 (e) representing an ablation of the cornea in the 

20 finite element analysis model by assigning values of 

21 thicknesses corresponding to the ablated portion of the cornea 

22 to nodes of the finite element analysis model in the ablated 

23 portion, respectively; 

24 (f) using the finite element analysis model, 

25 computing new values of the x,y,z coordinates at each of the 

2 6 nodes to simulate deformation of the cornea resulting from the 

27 proposed ablation; and 

28 (g) operating the computer to display the computed 

29 x,y,z coordinates to show the simulated deformation of the 

3 0 cornea . 

1 17 . A method of determining the modulus of elasticity of 

2 a surface of an eye at a plurality of locations of the 

3 surface, comprising the steps of: 

4 (a) storing in a computer system a finite element 

5 analysis model of an eye, the finite element analysis model 

6 including a number of nodes, the connectivities of which 

7 define a plurality of elements ; 

8 (b) applying a known external pressure to the eye 

9 and then measuring the topography of a portion of the eye 

10 using a topography measuring device to produce coordinates for 

11 a number of data points of the surface of the eye; 

12 (c) operating the computer system to interpolate 
• 13 between and extrapolate beyond the data points to obtain a 

14 reduced number of x,y,z coordinates that correspond to nodes 

' 15 of the finite element analysis model, respectively, and 

16 assigning the x,y,2 coordinates to the various nodes, 

17 respectively, and assigning the value of the external pressure 

18 to elements of the finite element analysis model corresponding 



Wp 94/18636 



PCT/US94/00908 



34 

19 to locations of the eye to which the known external pressure 

20 is applied in step (b) ; 

21 (d) determining a value representing intraocular 

22 pressure in the eye and assigning that value to each element; 

23 (e) assigning initial values of the modulus of 

24 elasticity to each element, respectively, of the finite 

25 element analysis model; 

26 (f) using the finite element analysis model, 

27 computing new values of the x,y,z coordinates at each of the 

28 nodes to simulate deformation of the cornea resulting from the 

29 external pressure and the intraocular pressure for the initial 

30 values of the modulus of elasticity; 

31 (g) operating the computer system to compare the 

32 new values of the x,y,z coordinates computed in step (f) with 

33 the x,y,z coordinates recited in step (c) ; 

34 (h) operating the computer system to modify values 

35 of the modulus of elasticity of the various elements of the 

36 finite element analysis model, respectively, if the comparing 

37 of step (g) indicates that a difference between the x,y,z 

38 coordinates of step (c) and the x,y,z coordinates computed in 

39 step (f) exceeds a predetermined criteria; and 

40 (i) repeating steps (f) through (h) until final 

41 values of the modulus of elasticity at each node are obtained. 



1 18. A method of determining change of a cornea of an eye 

2 as a result of an ablation on the cornea, the method including 

3 the steps of: 

4 (a) storing in a computer system a finite element 

5 analysis model of an eye, the finite element analysis model 

6 including a number of nodes, the connectivities of which 

7 define a plurality of elements; 

8 (b) applying a known external pressure to the eye 

9 and then measuring the topography of a portion of the eye 

10 using a topography measuring device to produce coordinates for 

11 a number of data points of the surface of the eye; 
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X2 (c) operating the computer system to interpolate 

13 between and extrapolate beyond the data points to obtain a 

. 14 reduced number of x,y,z coordinates that correspond to the 

15 nodes of the finite element analysis model, respectively, and 

16 assigning the x,y,z coordinates to the various nodes 

17 respectively, and assigning the value of the external pressure 

18 to elements of the finite element analysis model corresponding 

19 to locations of the eye to which the external pressure is 

20 applied in step (b) ; 

21 (d) determining a value representing intraocular 

22 pressure in the eye and assigning that value to each element; 

23 (e) assigning initial values of the modulus of 

24 elasticity to each element, respectively, of the finite 

25 element analysis model; 

26 (f) using the finite element analysis model, 

27 computing new values of the x,y,z coordinates at each of the 

28 nodes to simulate deformation of the cornea resulting from the 

29 external pressure and the intraocular pressure for the initial 

30 values of the modulus of elasticity; 

31 (g) comparing the new values of the x,y,z 

32 coordinates computed in step (f) with the x,y,z coordinates 

33 recited in step (c) ; 

34 (h) operating the computer system to modify values 

35 of the modulus of elasticity of the elements of the finite 

36 element analysis model, respectively, if the comparing of step 

37 (g) indicates a difference between the x,y,z coordinates 

38 obtained in step (c) and the x,y,z coordinates computed in 

39 step (f) exceeds a predetermined criteria; 

40 (i) repeating steps (f) through (h) until final 

41 values of the modulus of elasticity are obtained; 

42 (j) representing the ablation of the cornea in the 
• 43 finite element analysis model by assigning values of 

44 thicknesses corresponding to the ablated portion of the cornea 

' 45 to elements of the finite element analysis model in the 

46 ablated portion of the cornea, respectively; 

47 (k) using the finite element analysis model, 

48 computing new values of the x,y,z coordinates at each of the 
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49 nodes to simulate deformation of the cornea resulting from the 

50 proposed ablation; 

51 (1) comparing the simulated deformation of the 

52 cornea with a pre-established vision objective for the eye to 

53 determine if the ablation results in the vision objective 

54 being met; and 

55 (m) if the vision objective is not met, modifying 

56 the proposed ablation in the finite element analysis model and 

57 repeating steps (j) through (1) until the vision objective is 

58 met. 

1 19. A method of simulating change of a cornea of an eye 

2 as a result of a proposed incision on the cornea, the method 

3 including the steps of: 

4 (a) storing in a computer system a finite element 

5 analysis model of an eye, the finite element analysis model 

6 including a number of nodes, the connectivities of which 

7 define a plurality of elements; 

8 (b) applying a known external pressure to the eye 

9 and then measuring the topography of a portion of the eye 

10 using a topography measuring device to produce coordinates for 

11 a number of data points of the surface of the eye; 

12 (c) operating the computer system to interpolate 

13 between and extrapolate beyond the data points to obtain a 

14 reduced number of x,y,z coordinates that correspond to the 

15 nodes of the finite element analysis model, respectively, and 

16 assigning the x,y,z coordinates to the various nodes 

17 respectively, and assigning the value of the external pressure 

18 to elements of the finite element analysis model corresponding 

19 to locations of the eye to which the external pressure is 
2 0 applied in step (b) ; 

21 (d) determining a value representing intraocular 

22 pressure in the eye and assigning that value to each element; 

23 (e) assigning initial values of the modulus of 

24 elasticity to each element, respectively, of the finite 

25 element analysis model; 



wo 94/18636 



PCTAJS94/00908 



37 

2 6 (f) using the finite element analysis model, 
27 computing new values of the x,y,z coordinates at each of the 

. 28 nodes to simulate deformation of the cornea resulting from the 

29 external pressure and the intraocular pressure for the initial 

,30 values of the modulus of elasticity; 

31 (g) comparing the new values of the x,y,z 

32 coordinates computed in step (f) with the x,y,z coordinates 

3 3 recited in step (c) ; 

34 (h) operating the computer system to modify values 

3 5 of the modulus of elasticity of the elements of the finite 

36 element analysis model, respectively, if the comparing of step 

37 (g) indicates a difference between the x,y,z coordinates 

38 obtained in step (c) and the x,y,z coordinates computed in 

39 step (f) exceeds a predetermined criteria; 

40 (i) repeating steps (f) through (h) until a final 

41 value of the modulus of elasticity is obtained; 

42 (j) representing the incision in the finite element 

43 analysis model by means of a plurality of duplicate nodes 

44 adjacent to a group of incision-bounding nodes of the finite 

45 element analysis model located along the incision and 

46 representing the depth of the incision by means of a plurality 

47 of nonlinear spring elements each connecting an incision- 

48 bounding node to an adjacent duplicate node, respectively; 

49 (k) using the finite element analysis model, 

50 computing 

51 new values of the x,y,z coordinates at each of the nodes to 

52 simulate deformation of the cornea resulting from the incision 

53 and the intraocular pressure; 

54 (1) comparing the simulated deformation of the 

55 cornea with a pre-established vision objective for the eye to 

56 determine if the incision results in the vision objective 
'57 being met; and 

58 (m) if the vision objective is not met, modifying 

' 59 the incision in the finite element analysis model and 

60 repeating steps (j) through (1) until the vision objective is 

61 met. 
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1 20. A system for simulating deformation of a cornea of 

2 an eye in response to a corneal incision, comprising in 

3 combination: 

4 (a) a computer containing a finite element analysis 

5 model of an eye, the finite element analysis model including 

6 a predetermined number of nodes, the connectivities of which 

7 define a plurality of elements; 

8 (b) corneal topographer means for measuring the 

9 topography of a portion of the eye to produce coordinates for 

10 a number of data points of the surface of the eye; 

11 (c) means in the computer for interpolating between 

12 and extrapolating beyond the data points to obtain a reduced 

13 number of x,y,z coordinates that correspond to nodes of the 

14 finite element analysis model, respectively, and assigning the 

15 x,y,z coordinates to the various nodes, respectively; 

16 (d) means for measuring a value representing 

17 intraocular pressure in the eye; 

18 (e) means for inputting the value into the finite 

19 element analysis model and assigning that value to each 

20 element; 

21 (f) means in the computer for representing a first 

22 incision in the finite element analysis model by means of a 

23 plurality of duplicate nodes adjacent to a group of incision- 
• 24 bounding nodes of the finite element analysis model located 

25 along the first incision and representing a depth of the first 

26 incision by means of a plurality of nonlinear spring elements 

27 each connecting an incision-bounding node to an adjacent 

28 duplicate node, respectively; 

29 (g) means in the computer for executing a program 

30 including the finite element analysis model to compute new 

31 values of the x,y,z coordinates resulting from the first 

32 incision at each of the nodes, respectively; and 

33 (h) means responsive to the computer for displaying 

34 the nodes having the computed x,y,z coordinates to show the 

35 deformation of the cornea. 
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1 21. The system of Claim 20 including means in the 

2 computer for comparing the simulated deformation of the cornea 
. 3 with a pre-established vision oojective for the eye to 

4 determine if the first incision results in the vision 

5 objective being met, and means for modifying the first 

6 incision and/or adding a second incision in the finite element 

7 analysis model similar to the first incision if the vision 

8 objective is not met. 

1 22. The system of Claim 21 including means for measuring 

2 the thickness of various points of the cornea and/or sclera, 

3 and means for inputting the thicknesses into the computer and 

4 assigning values of the measured thicknesses to each element 

5 of the finite element analysis model, respectively. 

1 23. The system of Claim 22 wherein the modifying means 

2 includes a nonlinear programming computer program to determine 

3 how much to modify the number of incisions, the orientations 

4 of the incision, and the lengths and depths of the various 

5 incisions. 

1 24. A method of simulating deformation of a cornea of an 

2 eye in response to corneal incisions, comprising the steps of: 

3 (a) measuring the topography of a portion of the 

4 eye using a topography measuring device to produce coordinates 

5 for a large number of data points of the surface of the eye; 

6 (b) storing in a computer system a mathematical 

7 analysis model of the eye, the model including a number of 
^ 8 nodes, the connectivities of which define a plurality of 

9 elements ; 

'10 (c) determining a value representing intraocular 

11 pressure in the eye and assigning that value to each element; 

12 (d) representing the incision in the mathematical 

13 analysis model by means of a plurality of duplicate nodes 
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14 adjacent to a group of incision-bounding nodes of the 

15 mathematical analysis model located along the incision and 

16 representing the depth of the incision by means of a plurality 

17 of nonlinear spring elements each connecting an incision- 

18 bounding node to an adjacent duplicate node, respectively; and 

19 (e) using the mathematical analysis model to 

20 compute new values of the x,y,z coordinates resulting from the 

21 incision at each of the nodes, respectively. 

1 25. The method of Claim 24 including establishing a 

2 vision objective for the eye, wherein step (e) includes 

3 comparing the simulated deformation of the cornea with the 

4 vision objective to determine if the incision results in the 

5 vision objective being met, and, if the vision objective is 

6 not met, modifying the incision and/or adding another incision 

7 in the mathematical analysis model and repeating step (e) to 

8 determine if the vision objective is met. 

1 26. A method of simulating deformation of a cornea of an 

2 eye in response to a corneal incision, comprising the steps 

3 of: 

4 (a) measuring the topography of a portion of the 

5 eye using a topography measuring device to produce coordinates 

6 for a number of data points of the surface of the eye; 

7 (b) storing in a computer system a mathematical 

8 analysis model of the eye, the model including a number of 

9 nodes, the connectivities of which define a plurality of 

10 elements; 

11 (c) determining a value representing intraocular 

12 pressure in the eye and assigning that value to each element; 

13 (d) representing incision in the mathematical 

14 analysis model by assigning reduced thickness values and/or a 

15 reduced value of the modulus of elasticity to elements 

16 surrounding the incision; and 
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17 (e) using the mathematical analysis model to 

18 compute new values of the x,y,z coordinates resulting from the 
,19 incision at each of the nodes, respectively. 

1 27. A method of simulating deformation of a cornea of an 

2 eye as a result of a corneal ablation, comprising the steps 

3 of: 

4 (a) measuring the topography of a portion of the 

5 eye using a topography measuring device to produce coordinates 

6 for a number of data points of the surface of the eye; 

7 (b) storing in a computer system a mathematical 

8 analysis model of the eye, the model including a predetermined 

9 number of nodes, the connectivities of which define a 

10 plurality of elements; 

11 (c) determining a value representing intraocular 

12 pressure in the eye and assigning that value to each element; 

13 (d) representing the ablation of the cornea in the 

14 mathematical analysis model by assigning reduced values of the 

15 thicknesses and/or a reduced value of the modulus of 

16 elasticity to elements corresponding to the ablated portion of 

17 the cornea; and 

18 (e) using the mathematical analysis model to 

19 compute new values of the x,y,z coordinates resulting from the 

20 ablation at each of the nodes, respectively. 

1 28. The method of Claim 27 including establishing a 

2 vision objective for the eye, wherein step (e) includes 

3 comparing the simulated deformation of the cornea with the 

4 vision objective to determine if the ablation results in the 
" 5 vision objective being met, and, if the vision objective is 

6 not met, modifying the ablation in the mathematical analysis 

' 7 model and repeating step (e) to determine if the vision 

8 objective is met. 
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1 29. A method of simulating change of a cornea as a 

2 result of a proposed ablation on the cornea, the method 

3 including the steps of: 

4 (a) storing in a computer system a mathematical 

5 analysis model of an eye, the finite element analysis model 

6 including a number of nodes, the connectivities of which 

7 define a plurality of elements; 

8 (b) applying a known external pressure to the eye 

9 and then measuring the topography of a portion of the eye 

10 using a topography measuring device to produce coordinates for 

11 a number of data points of the surface of the eye; 

12 (c) operating the computer system to interpolate 

13 between and extrapolate beyond the data points to obtain a 

14 reduced number of x,y,z coordinates that correspond to the 

15 nodes of the finite element analysis model, respectively, and 

16 assigning the x,y,z coordinates to the various nodes 

17 respectively, and assigning the value of the external pressure 

18 to elements of the finite element analysis model corresponding 

19 to locations of the eye to which the external pressure is 

20 applied in step (b) ; 

21 (d) determining a value representing intraocular 

22 pressure in the eye and assigning that value to each element; 

23 (e) assigning initial values of the modulus of 

24 elasticity to each element, respectively, of the finite 

25 element analysis model; 

26 (f) using the mathematical analysis model, 

27 computing new values of the x,y,z coordinates at each of the 

28 nodes to simulate deformation of the cornea resulting from the 

29 external pressure and the intraocular pressure for the initial 

30 values of the modulus of elasticity; 

31 (g) comparing the new values of the x,y,2 

32 coordinates computed in step (f) with the x,y,z coordinates 

33 recited in step (c) ; 

34 (h) operating the computer system to modify values 

35 of the modulus of elasticity of the elements of the finite 

36 element analysis model, respectively, if the comparing of step 

37 (g) indicates a difference between the x,y,z coordinates 
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38 Obtained in step (c) and the x,y,z coordinates computed in 

39 step (f) exceeds a predetermined criteria; 

.40 (i) repeating steps (f) through (h) until a final 

41 value of the modulus of elasticity is obtained; 

42 (j) representing a proposed ablation of the cornea 

43 in the mathematical analysis model by assigning reduced values 

44 of thicknesses and/or a reduced value of the modulus of 

45 elasticity to elements corresponding to the ablated portion of 

46 the cornea; 

47 (k) using the finite element analysis model, 

48 computing new values of the x,y,z coordinates at each of the 

49 nodes to simulate deformation of the cornea resulting from the 

50 ablation and the intraocular pressure; 

51 (1) comparing the simulated deformation of the 

52 cornea with a pre-established vision objective for the eye to 

53 determine if the ablation results in the vision objective 

54 being met; and 

55 (m) if the vision objective is not met, modifying 

56 the ablation in the finite element analysis model and 

57 repeating steps (j) through (1) until the vision objective is 

58 met. 
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